Radiation efficiencies of the pulsars detected in the optical range 
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Using available multiwavelength data for the pulsars of different ages detected in the optical 
range, we analyze the efficiencies of the conversion of the pulsar spindown power E into the 
observed non-thermal luminosity L in different spectral domains. We find that these pulsars show 
significantly non-monotonic dependencies of the optical and X-ray efficiencies (77 = L/E) versus 
pulsar age with a pronounced minimum at the beginning of the middle-age epoch (~ 10^ yr) 
and comparably higher efficiencies of younger and older pulsars. In addition, we find a strong 
correlation between the optical and 2-10 keV X-ray luminosities of these pulsars that implies 
the same origin of their nonthermal emission in both spectral domains. 



1. Introduction 

The properties of the observed pulsar 
emission show that radio and gamma-ray 
photons have nonthermal origin and are 
produced by relativistic particles in dif- 
ferent parts of the pulsar magnetosphere. 
The soft X-ray-optical radiation can be a 
combination of several nonthermal spec- 
tral components from the magnetosphere 
and thermal components from the whole 
surface of a cooling neutron star (NS) 
and from its hot polar caps heated by 
the relativistic particles from the magne- 
tosphere. The nonthermal component has 
been found to dominate the optical emis- 
sion of several young, middle-aged, and 
old pulsars. It is believed that nonther- 
mal components are powered by NS rota- 
tional energy loss E = An^IP/P^, where 



/ 1.1 X 10^^ gm CIV? is the inertia mo- 
mentum of a NS of 1.4 Mq and 10 km 
radius, P and P are the pulsar period and 
its derivative, respectively. E is frequently 
called also the spindown luminosity and 
denoted as L^^. The ratio rj = L/E, where 
L is the photon luminosity, is used to char- 
acterize the efficiency of the conversion of 
E = Lsd into the pulsar emission at a 
given spectral domain. In this paper, we 
analyze the evolution of the nonthermal 
luminosities and efficiencies of pulsars de- 
tected in the optical range. Such analy- 
sis has been performed about 10 yr ago 
in pioneering works |ll2j . However, con- 
siderable progress in the pulsar distance 
measurements, and much higher quality 
of recent optical and X-ray observations 
allowing to resolve better the nonthermal 
and thermal emission components from 
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Table 1 

The dynamical ages r = P/2P, distances d, spindown luminosities E (or Lsd), and ob- 
served nonthermal luminosities in the radio, Lr, optical, Lopt, X-rays, Lx, and gamma- 
rays, L^, of seven radio pulsars detected in the optical range |4f 5j . Figures in brackets are 
±1(7 uncertainties of the values. 



Source 


log r 


d 


log E 


log Lr" 


log Lopt 


log Lx 


log 




[yr] 


[pc] 


[erg s^"*"] 


[mJy kpc^] 
408 MHz 


[erg s^-^] 
B band 


[erg s""*"] 
2-10 keV 


[erg s~"^] 
>400MeV 


Crab 


3.1 


2.0x10^ 


38.65 


3.41(4) 


33.23(5) 


36.67(1^°) 


35.7(1^) 


B0540-69 


3.2 


5.0x10^ 


38.17 


3.30(9) 


33.47(15) 


36.99(1^^) 


< 35.97 


Vela 


4.1 


293(11?) 


36.84 


2.64(20) 


28.3(3) 


31.2(11)^ 


33.9(1^) 


B0656+14 


5.0 


288(1^) 


34.58 


-0.27(9) 


27.53(8) 


30.30(11) 


32.37(1^°) 


Geminga 


5.5 


153(1^D 


34.51 


0.375(1^) 


26.95(11^) 


29.35(11) 


32.95(1^°) 


B1929+10 


6.5 


3611% 


33.59 


1.52(5) 


27.26(1^°) 


29.86(1^^) 


< 32.57 


B0950+08 


7.2 


262(5) 


32.75 


1.44(16) 


26.88(8) 


29.28(li^) 


< 32.51 



° The radio luminosity is defined traditionally as L^^ — S408d^ [mJy kpc^], where S408 is the observed flux 
density from a pulsar at 408 MHz in mJy, and di is its distance in kpc. At a typical radio band FWHM 
of about 100 kHz the conversion factor to standard luminosity units is « 9.51 x 10^^ [erg s~^]/[mJy kpc^]. 



the pulsars and to identify in these ranges 
several new objects enable us to update 
significantly previous analysis and to get 
qualitatively new results on the evolution 
of old pulsars. 

2. The sample and data sources 

We have analyzed multiwavelength non- 
thermal luminosities of seven pulsars de- 
tected in the optical range (see Table 1; 
|4|5j ). The advantage of this sample is a 
small (<10%) uncertainty in the distances 
to the pulsars, determined from the radio 
or optical parallax measurements. This 
strongly minimizes errors in the inferred 
luminosities. Distances for a few other 



optical pulsars are much less certain and 
we excluded them from consideration. In 
the optical we used the data obtained in 
the B-band |4p5j since practically all the 
pulsars from our sample were detected in 
this band and fluxes in other bands are 
not strongly different from these. Using 
this band also allows us to escape possi- 
ble contamination of the pulsar fluxes by 
strong nebular lines in cases when an as- 
sociated nebula around a pulsar exists, or 
is expected to exit. In X-rays we used 
the data obtained in 2-10 keV range |6|5j . 
where the nonthermal power law spectral 
tails dominate over the thermal compo- 
nents detected from several pulsars in a 
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Figure 1. From top to bottom: evolution 
of the pulsar spindown, radio, optical, X- 
ray, and gamma-ray luminosities and re- 
spective efficiencies, with dynamical age. 
Filled circles are the data taken from [7j. 
Hereafter errorbars and arrows show itlcr 
value uncertainties and la upper limits, 
respectively. 

softer energy range. We added also the 
data from T] on 5 old pulsars recently ob- 
served with the Chandra and XMM. In 
the radio and gamma-rays we used the 
data obtained at 400 MHz (ATNF catalog 
ID) and in the EGRET range of 400 MeV- 
10 GeV P, respectively. Here we con- 
sider only isotropic equivalents of the lu- 
minosities neglecting the emission beam- 
ing which is not yet properly known in the 
all considered ranges for all our objects. 




from Zavlin&Pavlov (2004) 
■ from PossenlL et al. (2002) without MSPs 
— for PSRs witli log T > 5.5 




Figure 2. Top: radio efficiencies of 500 
pulsars from ATNF [S] vs dynamical age; 
triangles mark the pulsars with more ac- 
curate parallax-based distances and circles 
select optical pulsars from Table 1. Bot- 
tom: X-ray efficiencies vs dynamical age; 
lines are the linear regression fits of evo- 
lution tracks of younger and older pulsars 
and numbers show the line slopes. 



3. Results 

3.1. Evolution 

The nonthermal luminosity and effi- 
ciency evolution of pulsars in different 
spectral domains, based on the selected 
sample of the optical pulsars, is shown in 
Fig. 1. We note the following points: 
(1) the spindown luminosity decreases 
monotonically with age, as expected from 
a formal dependence E ~ r^^P^^ and 



4 




Log r|_ I (B-band) 



Figure 3. Relationship between the opti- 
cal and X-ray efficiencies in the B band 
and 2-10 keV ranges. 
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Figure 4. The evolution of the optical and 
X-ray luminosities ratio with dynamical 
age. 

longer periods of older pulsars; (2) the 
optical, X-ray and gamma-ray luminosi- 
ties also decrease with age, but become 
almost constant starting from a middle- 
age epoch of 10'^ — 10^ yr; (3) the ra- 
dio luminosity is less monotonic but starts 
to increase gradually approximately from 
the same age; (4) result, the evolu- 
tion of the optical and X-ray efficiencies is 
non-monotonic and shows a pronounced 
minimum efficiency at a middle-age epoch 
of 10^ — 10'' years and comparable higher 
efficiencies for younger and older pulsars 
in contrast to the conception of a mono- 
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tonic decrease of the efficiency with age 
suggested by through lack of data on 
old pulsars at that time; (5) the shapes 
of the dependencies rjopt{T) and tix^t) are 
remarkably similar to each other; (6) the 
radio and gamma-ray efficiencies increase 
with age but their increases becomes ap- 
parently steeper from the same middle-age 
epoch. 

The fact that older pulsars are much 
more efficient radio emitters than younger 
ones is confirmed at a much higher signifi- 
cance level based on a rich sample of radio 
pulsars extracted from the ATNF catalog 
(Fig. 2, top panel). The scatter of the data 
along the mean evolution track shown in 
this panel by solid line may be explained 
by different geometrical/beaming factors 
and uncertainties in the distance. The 
pulsars with parallax based distances are 
generally below this line. This may re- 
flect a systematic shift between the paral- 
lax and dispersion measure distances used 
for the rest objects. 

The nonmonotonic efficiency evolution 
in X-rays becomes much less obvious if 
we consider a larger sample of pulsars 
identified in X-rays including a few of 
them identified in the optical and X-rays 
(Fig. 2. low panel). The statistics is much 
poor than in the radio and this picture is 
more affected by the distance and flux un- 
certainties. Nevertheless, applying linear 
regression ffis separately to sub-samples 
of sub-middle-aged and post-middle-aged 
pulsars, one can resolve hints of the non- 
monotonic evolution clearly seen in the 
optical sub-sample. In the example shown 
in Fig. 2, (low panel) the Geminga age of 
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IQ5.5 yj,^ when the optical efficiency start 
to increase, was used to divide by sub- 
samples. In this case the evolution track 
of older pulsars may be compatible with 
what is seen in the radio range (cf. dashed 
and solid lines) . Using other dividing ages 
from the range of the middle-aged epoch 
of 10^ — 10^ yr does not change the picture 
qualitatively though slopes values can be 
different, but they are always negative and 
positive for younger and older pulsars, re- 
spectively. In all of the X-ray plots we 
added X-ray data for 5 old pulsars de- 
tected recently with Chandra and XMM 
These data are likely to be consistent 
with the proposed evolution picture. To 
conffim or reject it better quality X-ray 
data and distance information are needed. 

3.2. Correlation between the optical 
and X-ray emission 

Similar shapes of the optical and X-ray 
evolution tracks suggest a correlation be- 
tween the luminosities in these domains. 
This is confirmed by the linear regression 
fit of the rjopt vs rjx distribution shown in 
Fig. 121 The correlation coefficient is 0.97. 
A strong correlation implies the same ori- 
gin for the emission in both spectral do- 
mains. It is independent of the differ- 
ence in the spectral slopes in these ranges 
observed for the pulsars of different age. 
For instance, the extrapolation of the non- 
thermal component of the X-ray spectrum 
of young Crab-like pulsars strongly over- 
shoots the observed optical fluxes, while 
for older pulsars the optical fluxes are gen- 
erally compatible with such extrapolation. 
This is reflected in the evolution of the ra- 
tio Lopt/Lx shown in Fig. HI The ratio 
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Figure 5. From top to bottom: relation- 
ships between the efficiencies in the radio, 
optical. X-ray, and 7-ray spectral domains 
and spindown luminosity. Different pul- 
sars are denoted by the same symbols, as 
in Fig. 1. Dashed and solid lines show the 
best linear regression ffis of the whole set 
of pulsars, and without the two younger 
ones (Crab and PSR B0540-69, with the 
highest Lsd), respectively. The slope val- 
ues are indicated near the lines. Dotted 
lines indicate a slope in case of the linear 
proportionality of the luminosity to the 
Goldreich- Julian current [Sj. 

becomes almost independent of age start- 
ing from a middle-age epoch. 

3.3. The efficiency and the spindown 
luminosity 

The distributions of the efficiencies of 
the optical pulsars in different spectral do- 
mains vs spindown luminosity are shown 
in Fig. El These distributions demonstrate 
that: (1) in the radio and gamma-rays. 
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Figure 6. Top: radio efficiencies of 500 
pulsars from ATNF database |H] vs L^^; 
notations are the same as in Fig. 2. Bot- 
tom: X-ray efficiencies vs spin-down lumi- 
nosity; lines fit the distribution of lower 
and higher Lgd pulsars; numbers show the 
line slopes. 

less energetic pulsars (with low Lsd) are 
more efficient photon emitters than more 
energetic ones; (2) in these ranges the ef- 
ficiency increases monotonically with the 
decrease; (3) the same is seen in the 
optical and X-ray ranges with an exclu- 
sion of the two most energetic and young 
pulsars, Crab and PSR B0540-69, which 
form a separate sub-class in these distri- 
butions; (4) recent data for the five old 
pulsars identified in X-rays (filled circles 
in Fig. 5) do not change the above ten- 
dency in X-rays and may even strengthen 
it (cf. Fig. ini bottom panel). 
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As in case of the evolution picture, 
the fact that the pulsars with lower Lsd 
are more effective radio emitters is ffimly 
confirmed by the analysis of a much 
larger sample of pulsars extracted from 
the ATNF catalog (Fig. El top panel). 
Comparing the triangles and black point 
distributions in Fig. IHl (top panel), we 
can note that the scatter around an av- 
erage distribution shown by the solid line 
is likely equally caused by uncertainties in 
the pulsar geometry (triangle points), and 
DM-based distances (black points). The 
slope of the distribution is different from 
a trivial one assuming ~ const (dashed 
line in this plot) but suggests an increase 
of the radio luminosity oc L^'J. 

The separation into two sub-classes of 
the low and high spindown luminosity pul- 
sars that is clearly seen in the middle 
panels of Fig. |31 for the optical pulsars 
is not obvious for a wider sample of pul- 
sars identified in X-rays (Fig. bottom 
panel). The reasons are likely to be the 
same as we have mentioned in Sect. 3.2 
(poor statistics + distance and flux un- 
certainties). However, as in the case of 
the evolution picture for the same sample 
shown in Fig. |3[ linear regression fits allow 
us to resolve a hint of a bimodal distribu- 
tion with the efficiency increases toward 
the low and high sides of the considered 
Lsd range. The bottom panel of Fig. 
where Lsd of Geminga was used, as in case 
of Fig. 121 to divide the pulsars by sub- 
samples, suggests also that the distribu- 
tion of the low Lsd pulsars may be roughly 
compatible with what is seen in the radio 
range (cf. solid and dashed lines in this 
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plot). The suggested bimodal distribution 
can hardly be explained by geometry ef- 
fects, if we compare the relatively small 
scattering of the triangle points around a 
mean distribution in the radio (top panel) , 
likely caused by geometry effects, and the 
much larger scattering of all points in X- 
rays (bottom panel). 

4. Summary 

Our multiwavelength analysis of the 
nonthermal luminosity and efficiency dis- 
tributions vs dynamical age and spin- 
down power of the pulsars detected in the 
optical range suggests significantly non- 
monotonic evolution of their emission in 
the optical and X-ray domains and/or a 
bimodal distribution of the respective effi- 
ciencies vs spindown power. Less energetic 
old pulsars are much more efficient pho- 
ton emitters than the middle-aged ones 
and their efficiency is comparable with 
that for the much younger and energetic 
Crab-like objects. The bimodal struc- 
ture is not significant but still resolved 
in a larger sample of pulsars identified 
in X-rays, where it may be smoothed by 
large distance and flux uncertainties and 
needs further study. An apparent evo- 
lution track of old pulsars in the optical 
and X-rays is roughly compatible with 
the track seen in the radio range where 
the evolution is accompanied by a signif- 
icant increase of the radio-efficiency with 
NS age. If it is true, this implies an uni- 
versal engine that drives the non-thermal 
evolution of old pulsars by a unified way 
increasing their efficiency with age in the 
all spectral domains. 
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